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ABSTRACT 

The luminosity dependence of kinematical and isophotal properties of elliptical galaxies is investi- 
gated using numerical simulations of galaxy merging, combined with semi-analytical models of hier- 
archical structure formation. Mergers of spiral galaxies as the only formation mechanism for elliptical 
galaxies can neither reproduce the kinematical and photometric properties of very massive elliptical 
galaxies nor the change from rotationally flattened disky to anisotropic boxy systems with increas- 
ing luminosity. We present numerical simulations showing that binary mergers of early-type galaxies 
open an additional channel for the formation of anisotropic, slowly rotating and boxy ellipticals. In- 
cluding this channel in a semi-analytical model we can successfully reproduce the observed trend 
that more luminous giant ellipticals are more boxy and less flattened by rotation. This trend can 
be strengthened by suppressing residual gas infall and star formation for galaxies with stellar bulge 
masses M, > 3 x IO^^Mq. Hence we propose that mergers of early-type galaxies play an important 
role for the assembly of massive elliptical galaxies. 

Subject headings: galaxies: elliptical - galaxies: interaction- galaxies: structure - galaxies: evolution 
- methods: numerical 



1. INTRODUCTION 

Giant elliptical galaxies show a strong dependence of 
their kinematic al and photometric propertie s with lumi- 
nosity (see e.g. iKormendv fc Bendedll99(l) . Tradition- 
ally the kinematics supporting the shape of ellipticals 
is parametrized by the anisotropy parameter (wmaj/co)* 
which relates the observed ratio of the velocity along the 
major axis and the central velocity dispersion, t^maj/coi 
to the value expected for an isotropic rotator v/a = 
•\/e/(l — e) with the observed ellipticity e l)Binnevlll978j) . 
Observed slow rotating (umaj/tro < 0.2) ellipticals in 
general have (t'maj/fo)* < 0.5 (see FigQ and their 
shape is assumed to be supported by anisotropic veloc- 
ity dispersions. In the following we call these systems 
anisotropic (A). They tend to be the brightest ellipti- 
cals with boxy isophotes and have the oldest, most metal 
rich stellar populations. The bul k of their stars forme d 
on very short t i me-sc ales (see e.g. iHeavens et all ()2004j) : 
IThomas et all l|2005j) and references therein). Low lu- 
minosity ellipticals rotate faster (wmaj/fo > 0.2) and 
have (wmaj/fo)* > 0.5. They are assumed to be flat- 
tened by additional rotation (R), have disky (sometimes 
boxy) isophotes, younger stellar populations, and formed 
on longer time-scales. Recently, observations as well 
as simulations have shown that disky ellipticals show a 
significant under ly ing velocity anisotr opy (Binnev 2005; 
iBurkert fc Naabli200l iCappellari et'al. 1 1200,^11) . In this 
sense remnants with (i^maj/co)* > 0.5 might be partly 
supported by anisotropic dispersions and partly by in- 
trinsic rotation. 

What is the origin of the strong mass dependence of 
the properties of elliptical galaxies? Does it fit into a 
scenario where ellipticals result from major galaxy merg- 
ers? While the properties of low luminosity ellipticals 
are consistent with being remnants of early-type spiral 
mergers (Naab & Burkert 2003(NB03); Naab & Truiillij 
120051) it is for several reasons unlikely that high lumi- 



nosity ellipticals have formed in a similar way: (i) the 
stellar populations of spirals are too young and metal 
poo r and their formation timescales are too long (see 
e.g. iBell et all 120051: IHeavens et al.ll200J and references 
therein), (ii) very luminous ellipticals are significantly 
more massive than typical spirals, (iii) numerical simula- 
tions indicate that kinematic and photometric properties 
of spiral merger remnants disagree wit h the properties 
of th e most massive ellipticals (NB03; iNaab fc Truiilld 
l2005j) . Although equal-mass spiral-s piral mergers can 
form slowly rotating, boxy remnants ijNaab et alJll99!j|) 
they do, in general, not result in a homogenous family 
si milar to bright ellipti c als. 

Khochfa r fc BurkertI 1)2003 ?) and Khochfar & BurkerJ 
(2005) (hereafter KB05) used semi-analytical model- 
ing to demonstrate that massive ellipticals should have 
assembled predominantly in mergers between bulge- 
dominated gas poor (dry) early- type galaxies. They 
concluded that the observed characteristic properties of 
massive ellipticals might reflect the physics of early-type 
mergers which in this case should lead to predominantly 
boxy and slowly rotating, anisotropic systems, indepen- 
dent of the progenitor mass ratio. This would be in 
contrast to spiral-spiral mergers which produce boxy, 
anisotropic systems only if the progenitors have equal 
masses. The semi-analytical models of KB05 predict that 
early-type mergers also occur at low redshifts. These 
merger remnants would still appear old if gas infall and 
subsequent star formation was negligible during and after 
the merger events. 

Some numerical simulations agree with the con- 
clusions of KB05. Gas- free mergers of spheroids 
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remnants are in general similar to observed ellipticals 
ijNaab fc BurkertI l2000t iGonzalez-Garcia fc van Albada! 
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However, a detailed investigation of the isophotal 
shape and kinematics of dry, early-type merger remnants 
was missing up to now. 

There is direct observational evidence for the 
existence of e arly-type mergers in gala x y clusters 
I van Dokkum et al. 19991 iTran et alJ l2005() . iBell et all 
( 2005a(l and Ivan DokkumI ( 2005) also find clear evidence 
for this process in the field. Comparing with merger 
timescales taken from the simulations presented here 
IBell et all l|2005aj) conclude that a spheroidal galaxy 
should have undergone typically 0.5-2 early-type major 
mergers since redshift z k, 0.7. Investigations of the 
evolution of the luminosity function of elliptical galax- 
ies using various surveys show that the number of early- 
type galaxies with luminosities 0.5-2 has been growing 
since a redshift oi z k, \ accompanied by an increase of 
the total stellar mass in the early-ty pe galaxy population 
l)Bell et al.ll2"00l iFaber et al.ll2005l() . These observations 
support the scenario that massive galaxies assemble a 
significant amount of their mass even at z < 1 by ma- 
jor mergers of early- type galaxies and accretion of satel- 
lites rather than by s tar formation l)Drorv et alJ 1200^ 
IKhochfar fc "sTil3l200l . 

In this letter we present the kinematic and photometric 
properties of simulated early-type merger remnants as 
a function of their initial conditions. We combine the 
outcome of our simulations with semi-analytical models 
and predict the present day abundance of anisotropic and 
rotationally flattened ellipticals. 

2. THE MERGER MODELS 

As prototypical progenitors for mergers between early 
type galaxies we use an anisotropic remnant of a spiral- 
spiral merger with a mass ratio 1:1 and a rotationally flat- 
tened remnant of a 3:1 spiral-spiral merger from the sam- 
ple of NB03. These objects resemble low and intermedi- 
ate mass elliptical galaxies with respec t to their photo- 
metri c and kinematic properties fNB03: lNaab fc Truiillcl 
l2005^ and therefore can be considered as very good, but 
not perfect models for early-type galaxies. However, the 
mismatch in higher order moments of the linc-of-sight 
velocity distribution of 3:1 remnants and observed ellip- 
ticals (jNaab fc Burkert 2001 '! does not influence the con- 
clusions based on global kinematics of early-type merger 
remnants presented here. We performed second gener- 
ation binary mergers with mass ratios of = 1 and 
77 = 3 between anisotropic (llAA, 31AA) and rotating 
(IIRR, 31RR) progenitors. The low-mass progenitor in 
the unequal-mass mergers contained a fraction of I/77 
of the mass and particle number of the high-mass com- 
panion in each subcomponent and its size was scaled as 
yJTJrj. Particle numbers were the same as in NB03, re- 
sulting in 360000 luminous and 640000 dark particles for 
llAA. 

Test simulations indicated that the initial orientation 
of the early-type progenitors is less important than for 
spiral-spiral mergers. Therefore we decided to focus 
on a variation of the pericenter distances of the merg- 
ers as an analysis of cold dark matter simulations by 
IKhochfar fc Burkerti l)2005a|) has shown that impact pa- 
rameters of major mergers are larger than typically as- 
sumed in isolated simulations. For every initial setup 
the galaxies approached each other on nearly parabolic 
orbits with three pericenter distances of rgop = 0-5, 2, 6 



where the disk scale length of the more massive first 
generation spiral was hd — I (see NB03 for all details) 
resulting in 12 remnants in total. The initial separa- 
tion was fixed to 30 length units. The simulations were 
performed with the treecode VINE in combinatio n with 
the s pecial purpose hardware GRAPE-5 (K awai et alJ 
|2000() . The gravitational Plummer-softening length was 
the same (cpium — 0.05) as for the first generation spi- 
ral mergers. The particle orbits were integrated with a 
leapfrog integrator and a fixed time-step of At = 0.02 
which corresponds to « 0.2% of the half-mass rotation 
period of the first generation disks. The total energy 
was conserved better than 0.6%. Every remnant was 
analyzed 10 dynamical timescales after the merger was 
complete. 

The isophotal shape and the kinematics of 500 ran- 
dom projections of each of the 12 merger remnants were 
analyzed. For each projection the characteristic elliptic- 
ity eoff and isophotal shape parameter a4efF as well as 
the rotational velocity along the major axis, Umaj, and 
the central velocity dispersion, ctq, was determined as in 
NB03. 

The left panel of FigQ shows the location of the rem- 
nants in the (wmaj/fo)* — 04 plane. For every merger 
pair the projected remnants of the three orbital ge- 
ometries have been combined. All llAA (solid line), 
llRR (dashed line) and 31AA (dotted line) remnants 
have (wmaj/co)* < 0.5 and show predominantly boxy 
isophotes. Only 31AA remnants with a large pericen- 
ter distance, rgop = 6 and llRR remnants with a small 
pericenter distance rgop — 0.5 can appear significantly 
disky (a4off > 0.5). 31RR remnants show distinct prop- 
erties and have (wmaj/fo)* > 0.5 with predominantly 
disky isophotes. 

The distinct kinematical properties of the merger rem- 
nants becomes more evident in the (wmaj/co) — eoff-plane 
which is shown in the right panel of Fig^ llAA, 
llRR, and 31AA remnants are slow rotators in good 
agreement with observed boxy and anisotropic ellipticals 
(filled squares). 31RR remnants show stronger rotation 
(0.2 < fmaj/fo < 0.5) and are located in the overlap 
region between boxy and disky ellipticals. This analysis 
indicates that unequal mass mergers of early-type galax- 
ies can result in slowly rotating anisotropic systems if the 
more massive progenitor had similar properties. 

3. SIGNIFICANCE OF EARLY-TYPE MERGERS AND THE 
MASS DEPENDENCE OF ELLIPTICAL GALAXY 
PROPERTIES 

Semi-analytical models are useful to quantify the im- 
portance of early-type mergers for elliptical galaxy for- 
mation and to compare the theoretical models with ob- 
servations. The model used here is described in detail in 
KB05. For the present study we adopt the following cos- 
mological parameters: fig = 0.3, J^a = 0.7, h = 0.65 and 
fJs = 0.9. Early type galaxies were defined by a bulge to 
total mass fraction of > 60%. All galaxies with smaller 
bulge fractions were defined as spirals. As the last major 
merger determines the final properties of the elliptical 
galaxies we show in Fig|21the predicted number ratio of 
last elliptical-elliptical (E-E) mergers to the sum of the 
last elliptical-spiral (E-S) and the last spiral-spiral (S-S) 
mergers found in our simulations. The importance of 
E-E mergers steadily increases with increasing remnant 
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Fig. 1. — Left panel: Isophotal shape aA^s versus anisotropy parameter (fmaj /""O ) * ■ The contours indicate the 80% probability to find a 
remnant of 1:1 anisotropic (llAA; solid), 1:1 rotationally flattened (IIRR; dashed), 3:1 anisotropic (31AA; dotted), and 3:1 rotating (31RR; 
dashed-dotted) early-type mergers in the enclosed area. For every initial galaxy pair the remnants of the three adopted impact parameters 
have been combined. The horizontal line at (fmaj /""o)* = 0-5 separates the anisotropic boxy from the more rotationally flattened disky and 
boxy ellipticals. Right panel: EUipticity e^ff versus the ratio of major axis rotation velocity and central velocity dispersion (fmaj/fo)- For 
comparison, in both panels the location of 1:1 (bright grey) and 3:1 (dark grey) spiral merger remnants are indicated by the shaded areas 
(NB03). The data for observed disky (diamonds), very anisotropic boxy (filled boxes) and less anisotropic boxy (open boxes) ellipticals 
have been kindly provided by Ralf Bender. IIAA, llRR, and 31AA remnants are consistent with slowly rotating, anisotropic and boxy 
elliptical galaxies. 




11.0 
log (M./Me) 

Fig. 2. — Number ratio of last major mergers between two el- 
liptical galaxies (E-E) and the sum of mergers between two spirals 
(S-S) and a spiral and an elliptical (E-S) as a function of stellar 
mass. The solid line shows the global average while the dotted 
and long dashed lines show results field- and high density (cluster) 
environments, respectively. The symbols show the maximum pos- 
sible remnant mass formed by E-E mergers at different rcdshifts 
for which S-S and E-S merger still occur. 

mass. Remnants with stellar masses > 6 X 1O"M0 
are generally assembled by E-E mergers. Our simulations 
predict that the Ne^e/{Ne-s + Ns-s) ~ Af,— relation 
shown in Fig [21 is independent of redshift and environ- 
ment. The only difference is that ellipticals reach larger 
masses in denser environments and that at higher red- 
shifts the relation stops at lower remnant masses. 

A summary of the conditions, which are based on the 
numerical simulations presented here and in NB03, that 
lead to anisotropic systems or systems with additional 
rotation is shown in Tab. ?? where S denotes spiral 
galaxies. For example the combination 1R-3A denotes 
a merger of a rotationally flattened elliptical with an 
anisotropic elliptical that is 3 times more massive. This 
merger is assumed to lead to an anisotropic (A) early- 
type galaxy. All simulations so far show that equal-mass 
spiral-spiral or elliptical-elliptical mergers predominantly 
lead to slowly rotating anisotropic systems. It is there- 
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Fig. 3. — Ratio of anisotropic ellipticals to the total number of 
elliptical galaxies as a function of their _B— band luminosity 
Stars show observational data provided by Ralf Bender. The dot- 
dashed line shows the result for S-S mergers alone which can not re- 
produce the observed trend. The SAM predictions including early- 
type mergers without (dashed line) and with (dotted line) assuming 
a shut-off mass scale follow the otjserved trend. The shaded areas 
indicate the scatter in the SAMs. 

fore very likely that this will also be true for equal-mass 
mixed mergers between spirals and ellipticals. In the case 
of unequal-mass mixed mergers we make the conserva- 
tive assumption that the property of the more massive 
progenitor dominates the structure of the remnant. We 
now can classify all the merger remnants in our semi- 
analytical models using the approach of KB05 and Table 
??. The result (dashed line) is compared to observations 
kindly provided by Ralf Bender ^ in FigOl As anisotropic 
systems also form in some unequal-mass early-type merg- 
ers the observed trend with galaxy luminosity can be 
reproduced. An even better agreement (dotted Hue) is 
found if gas infall and star formation is suppressed above 
a critical bulge mass of M* > 3 x lO^^Mf?) , moti vated 
bv recent predictions of e.g.~lBun dv et aLl |)200.'tD and 
iBell et al.l (^2005). In this case, the growth of stellar 
disks by gas infall and star formation is suppressed above 

extended data set based on IBender et a!] IT993) 



4 



this critical mass limit, increasing the frequency of early- 
type mergers and by this enhancing the formation rate 
of anisotropic systems. A model which generalizes the 
S-S results of NB03 to all mergers (dot-dashed line) fails 
to reproduce the observed trend. 

4. CONCLUSIONS 

Binary major mergers of spiral galaxies alone do not 
provide a viable formation mechanism for the whole pop- 
ulation of giant elliptical galaxies. In particular they fail 
to reproduce the dominance of anisotropic, slowly ro- 
tating and boxy ellipticals at high luminosities. A re- 
vised major merger scenario which in addition allows 
major mergers of gas poor early-type galaxies helps to 
solve this problem. Individual numerical simulations of 
coUisionless binary early-type mergers presented in this 
letter reveal an additional channel for the formation of 
anisotropic, slowly rotating ellipticals. In the context of 
hierarchical cold dark matter cosmologies this additional 
formation mechanism can successfully reproduce the ob- 
served trend that more luminous giant ellipticals are 
slower rotators and more supported by anisotropic veloc- 
ity dispersions alone. A trend that matches observations 
nicely, if residual gas infall, star formation (e.g. the sub- 
sequent formation of disks) for galaxies with stellar bulge 
masses > 3 x lO^^M© is suppressed. A threshold 
value like this has been proposed recently l|Biindv et al.l 
l2005tlBell et"am2nol . We can also conclude that dissi- 
pation becomes less important for the assembly of more 



massiv e ellipticals as proposed by iKormendv fc BendeJ 

The simulations also indicate that early-type merger 
remnants do pref erentially have boxy isophotes (see also 
iGovernato et ah * 1993). However, individual merger ge- 
ometries also lead to disky and anisotropic remnants. 
Further investigations will show whether this higher or- 
der effect requires the inclusion of additional physics like 
e.g. the coalescence of massive black holes. 

Our analysis strongly supports a scenario where ma- 
jor mergers between ellipticals play an important role 
for the assembly of elliptical galaxies even at low red- 
shifts. Provided that the stellar population of the pro- 
genitor galaxies was old even the late assembly of an 
elliptical galaxy in a major early-type merger will not 
change this property. Further investigations will have to 
show in how far this scenario is in details consistent with 
the redshift evolution of ellipticals and local scaling rela- 
tions like the color-magnitude relation, the Fundamental 
Plane and the black-hole-mass-tr relation. 
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